Magnetic and Orbital Order in (i?Mn03)„/ (AMn03)2 n Superlattices Studied via a 

Double-Exchange Model with Strain 
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The two-orbital double-exchange model is employed for the study of the magnetic and orbital 
orders in (_RMn03)„/(AMn03)2n {R- rare earths; A: alkaline earths) superlattices. The A-type 
antiferromagnetic order is observed in a broad region of parameter space for the case of SrTiOa 
as substrate, in agreement with recent experiments and first-principles calculations using these 
superlattices. In addition, also a C-type antiferromagnetic state is predicted to be stabilized when 
using substrates like LaAlOa with smaller lattice constants than SrTi03, again in agreement with 
first principles results. The physical mechanism for the stabilization of the A- and C- magnetic 
transitions is driven by the orbital splitting of the x 2 — y 2 and 3z 2 — r 2 orbitals. This splitting is 
induced by the Q3 mode of Jahn- Teller distortions created by the strain induced by the substrates. 
In addition to the special example of (LaMn03) n /(SrMn03)2n, our phase diagrams can be valuable 
for the case where the superlattices are prepared employing narrow bandwidth manganites. In 
particular, several non- homogenous magnetic profiles are predicted to occur in narrow bandwidth 
superlattices, highlighting the importance of carrying out investigations in this mostly unexplored 
area of research. 

PACS numbers: 71.30.+h, 73.21. Cd,75.47.Lx 
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I. INTRODUCTION 

In recent years, remarkable progress has been achieved 
in the experimental and theoretical study of com- 
plex oxide heterostructures. At present, atomic- 
scale smooth interfaces involving pairs of oxides with 
similar crystal structures and lattice constants can 
be routinely synthesized and characterized with so- 
phisticated techniques, potentially leading to a vari- 
ety of interesting electronic devices^— As a proto- 
type of correlated electron system, manganites have 
been widely involved in complex oxides heterostruc- 
tures, e.g. manganite-cuprate^^ manganite-titanate 
manganite-nickelate, 11 and mang anite-BiFe0 3 .- 12 Even 
pure manganite-manganite superlattices made by anti- 
ferromagnetic insulators LaMnC>3 (LMO) and SrMnC>3 
(SMO) are quite nontriviah— ~— A metal-insulator tran- 
sition was observed with increasing n for the case 
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of (LaMn03)2n/(SrMnC>3) n superlattices 



17.18.20 



with 



LMO to SMO in proportion of 2 to 1. Both employing 
microscopic models and density-functional theory (DFT) 
calculations, the presence of magnetic modulations in 
these superlattices have been reported, which may be 
responsible for the metal-insulator transition! 25 ' 26 

Recently, May et at observed an enhanced Neel tem- 
perature of the A-type antiferromagnetic ( A- AFM) order 
in the (LaMn0 3 )„/(SrMn03)2„ (n=l, 2) superlattices^ 



where now the LMO and SMO are in proportion of 1 
to 2. This is surprising since from the theory perspec- 
tive it is not straightforward to understand the pres- 
ence of a uniform A-AFM order because previous the- 
oretical studies (using both model and DFT approaches) 
predicted that the ground state spin order within the 
SrMnOs region was G-type antiferromagnetic (G-AFM) 
even when there were only two consecutive SrO sheets 
as in the case of (LaMn0 3 )4/(SrMn0 3 )2-2 5 .>2 6 In the 
(LaMn03)„/(SrMn03)2n superlattices, the consecutive 
SrO sheets can reach up to 4 when n=2. Thus, the e g 
electron density deep within this SrMnOs region is ex- 
pected to be very close to zero as in the case of bulk 
SrMnOs. For such a low electronic density, typically the 
G-AFM order is robust and the A-AFM order would be 
unusual according to the corresponding phase diagrams 
for bulk compounds. 27 However, contrary to these expec- 
tations, the experiments of May et al. unveiled a robust 
A-AFM order which was uniform all through the super- 
lattices without any noticeable magnetic modulation. 21 
Thus, it is very important to perform better refined the- 
oretical studies to understand this puzzling behavior. 

Very recently, some of the authors performed a 
DFT calculation that actually confirmed the presence 
of a nearly uniform A-AFM phase for the case of 
(LaMn0 3 )„/(SrMn0 3 )2n when using SrTi0 3 (STO) as 
a substrate and employing an on-site Hubbard U repul- 
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sion not too large (e.g. U < 1.5 eV)^ Furthermore, 
a uniform C-type antiferromagnetic (C-AFM) order was 
predicted to occur for these superlattices for the case of 
a LaA10 3 (LAO) substrate also for U < 1.5 eV^ In 
this previous DFT study, the effect of strain was em- 
phasized to tune the orbital and magnetic orders. The 
role of strain to stabilize A-AFM and C-AFM phases 
for Lai_ x Sr x Mn0 3 thin him and (LaMn0 3 )i/(SrMn0 3 )i 
superlattice was also discussed before from the perspec- 
tive of DFT calculations^^ 

In the present manuscript, the (i?Mn0 3 ) n /(^4Mn0 3 )2n 
superlattices will be revisited now using the two-orbital 
double-exchange model, as opposed to a first-principles 
study. Our mission is not only to confirm (or refute) 
the DFT predictions from a conceptually different per- 
spective, but also to reveal in more detail the physical 
mechanisms behind the results, particularly with regards 
to the relationship between magnetism and orbital pop- 
ulation. Furthermore, our model calculations can go be- 
yond (LaMn0 3 )„/(SrMn0 3 )2n and shed light on generic 
narrow-bandwidth (i?Mn0 3 ),i/(AMn0 3 )2n superlattices 
which have not been experimentally prepared thus far 
and that are difficult to analyze accurately with DFT 
methods. Our study confirms the presence of the A-AFM 
and C-AFM states, in good agreement with DFT, but it 
has also unveiled a plethora of other phases, including 
exotic states with non-uniform configurations, that po- 
tentially can emerge in narrow bandwidth superlattices 
after the appropriate selection of elements in the chemical 
formula and of the substrate. 

The organization of the manuscript is as follows. In 
Section II, the model and the method used are presented, 
with a discussion of the profile of the electrostatic poten- 
tial used for the superlattices as well as the variational 
states employed. The importance of using twisted bound- 
ary conditions, particularly for the A-AFM and C-AFM 
phases, is discussed in this section as well. Section III 
contains the main results and specific conclusions for the 
cases of the two substrates STO and LAO under study 
here. Emphasis is given to phase diagrams varying the 
superexchange and the electrostatic potential strength. 
Also in this section, the orbital occupation and electronic 
density for the many layers are discussed, as well as the 
implications of our results for more narrow bandwidth 
manganites. A summary is provided in Section IV. 



II. MODEL, METHOD, AND BOUNDARY 
CONDITIONS 



In this investigation, the standard two-orbital double- 
exchange model will be employed, in the widely used limit 
of an infinite Hund coupling for simplicity. Work car- 
ried out over several decades by dozens of groups has 
repeatedly shown that this double-exchange model pro- 
vides a realistic starting approximation for manganite 
compoundsi^Ir— More explicitly, the model Hamiltonian 



is given by 
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In this Hamiltonian, the first term represents the 
standard double-exchange hopping process for the e g 
electrons. The operators Ci a (cj Q ) annihilate (create) 
an e g electron at the orbital a of the lattice site i. 
The spin index is no longer necessary since the spin of 
the e g electrons is always parallel to the localized t2 g 
spin Si due to the infinite Hund coupling approxima- 
tion. The three nearest-neighbor (NN) hopping direc- 
tions are denoted by r. The Berry phase f2y gener- 
ated by the infinite Hund coupling approximation equals 
cos(6>i/2) cos(6»j/2) + sin(6»,/2) sm{9j/2) exp[-i(<j) l - 
where 9 and </> are the polar and azimuthal angles of the 
t 2g spins, respectively, and i and j are NN sites. Two e g 
orbitals (a: d x i_ y 2 and b: d 3z 2_ r 2) are involved in the 
double-exchange process for manganites, with the hop- 
ping amplitudes given by: 
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Below, in the presentation of our results to will be con- 
sidered as the unit of energy. The second term of the 
Hamiltonian is the antiferromagnetic superexchange (SE) 
interaction between the NN t2 g spins. The typical value 
of the SE coupling Jafm is approximately between 0.05io 
and O.lto, based on a variety of previous investigations 
for bulk manganiteSi 25 i 31 ~ — The third term stands for the 
electron-lattice interaction, with A being a dimensionless 
coupling. A crucial aspect of the present investigation is 
the inclusion of the strain caused by the substrate. This 
strain is here modeled using the Q% (~ (25 z — 5 X — 5 y )) 
mode of the Jahn- Teller ( JT) distortions, where S a stands 
for the change in the length of the O-Mn-0 bonds along 
a particular axis a. The operator r 2 (— c^c a — c\cb) is an 
orbital pseudospin operator. The term AQ 3 r 2 splits the 
energy levels of the orbitals d x 2_ y 2 and d 3z 2_ r 2. 

From a formal perspective, the double-exchange hop- 
ping amplitudes, as well as the superexchange and Jahn- 
Teller couplings, may be different at the i?Mn0 3 and 
AMn0 3 regions, and at the interface. However, since 
the in-plane lattice constants are fixed to those of the 
substrate by the growing process, then the in-plane in- 
teraction couplings will be considered almost uniform all 
through the superlattices. With regards to the out-of- 
plane coupling strengths, in principle they should take 
different values for i?Mn0 3 , AMn0 3 , and the interface. 
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However, to fully consider these effects additional param- 
eters are needed since their precise values are not avail- 
able in the current state of investigations in oxide super- 
lattices. Thus, in order to avoid having too many free pa- 
rameters, it is a practical choice to carry out these inves- 
tigations assuming uniform couplings, as used in previous 
work. (25,34,35) This approximation reduces the numbers 
of model parameters needed in our study, allowing us to 
focus on the most important aspects of the physics of the 
problemj 25 i 34 ' 35 This approximation reduces the numbers 
of model parameters needed in our study, allowing us to 
focus on the most important aspects of the physics of the 
problem. 

In our study neither the breathing mode Qi nor the JT 
mode Q2 are considered in the calculations for the follow- 
ing reasons. First, the contribution that arises from Q\ 1 
which effectively modifies the on-site potential, can be 
simply merged into the fourth term. Second, with re- 
gards to the Q2 JT mode its average value < Q2 > is 
zero if the substrate is biaxially isotropic. Although the 
local values of Qn may become nonzero in special cases 
such as for the undoped manganites i?MnC>3 with or- 
bital order ^ the Qn profile is not considered important 
in the bulk doped manganites such as i?i/3A 2 /3Mn03 37 
and, more importantly for our studies on superlattices, 
they have not been reported as being of relevance in pre- 
vious studies of (i?Mn03)„/(AMn03)2n superlattices to 
our knowledge. Therefore, the Qn term will be neglected 
in the present effort for mere simplicity. 

In our calculations, Q-^ is proportional to the local 
quantity (c/a — 1), where c and a are the out-of-plane 
and in-plane lattice constants, respectively. The specific 
generic values of Q3; used in the present work are listed 
in Table |U In the present effort, the superlattices are as- 
sumed to be fully strained following the underlying sub- 
strates, as observed experimentally. 21 

The fourth term is the on-site Coulomb potential en- 
ergy, and rii is the local e g electronic density. Vi is the 
Coulomb potential arising from the A-site cations, from 
the e g electrons, and from the oxygen anions (i.e. the 
Qi mode distortion mentioned before). In the rest of the 
publication, the profile of Vj will be approximately deter- 
mined by the NN A-site cations ; 25 ' 38 as shown in Table U 
while its actual magnitude V will be considered a param- 
eter and phase diagrams will be constructed varying this 
parameter. Fixing such a profile is not a drastic approx- 
imation since previous more detailed studies, even us- 
ing the self-consistent Poisson equation, have shown that 
these sophisticated approaches tend to produce smooth 
potentials with profiles similar to those used here^ For 
the benefit of the reader, the distribution of the Mn02, 
RO, and AO layers present in our superlattices are shown 
in Fig. [U 

The above described model Hamiltonian has been nu- 
merically solved here using finite clusters. To simulate 
the case of the (RMn0 3 ) n /(AMn0 3 ) 2 n (n=l, 2) super- 
lattices, a 4 x 4 x 6 cluster with twisted boundary con- 
ditions (TBC)22 is here adopted. The TBC with proper 



■■■■ iiiii 1 

A Mn A Mn R Mn A Mn A Mn R Mn 

^ — ■ n=2 — ^= ^_ ^_ 





n 




2 


















































A Mn 




A 


Mn 


A 


Mn 




R 


Mn 




R 


Mn 


A 


Mn 

































1 2 T 3 . , 4 5 6 



FIG. 1. (Color online) Sketch showing the distribution of 
layers in the finite cluster with six Mn-oxide layers used in 
our studies. "Mn" stands for the Mn02 layers, "R" for the 
-RO layers, and "A" for the AO layers. The only active layers 
for the mobile electrons are the Mn layers in our study, but the 
other layers influence on the Mn layers by their electrostatic 
potential, as discussed in the text and shown in Table I. 



/c-meshes (kxkxk used here) are often employed to reach 
larger lattices, thus reducing the finite-size lattice effects 
(which are indeed of relevance in the present study, as 
shown below).— 

To determine the ground state phases in these superlat- 
tices, a variational procedure is used, namely the energies 
of candidate phases are compared in the Jafm-^ param- 
eter space, for several values of XQ3 as an additional pa- 
rameter. The candidate phases are proposed according 
to the symmetries of the superlattices. Contrary to other 
variational studies limited to a relatively small number 
of proposed states, in this effort 30 candidates for the 
n=l case and 56 candidates for the n—2 case have been 
considered. These candidate states were constructed by 
combining the well-known FM, A-AFM, C-AFM, and G- 
AFM spin arrangements in each of the layers (note that 
the CE-AFM state is not considered since the Q 2 mode is 
neglected). Some of these candidates will appear explic- 
itly in the phase diagrams described below in Section III, 
and some of them are presented in Fig. [2] for the benefit 
of the reader. 



A. Discussion of finite-size lattice effects 

In several previous computational studies of the 
double-exchange model, periodic boundary conditions 
(PBC) were used for Monte Carlo simulations carried 
out in real-space lattices.— Due to the substantial CPU- 
time consumed by the frequent diagonalizations in the 
fermionic sector, calculations have been typically lim- 
ited to relatively small lattice sizes, especially for three- 
dimensional studies. For example, in several cases the in- 
plane lattice cell was only 4 x 4 in size ) 25 ' 35 i 41 which might 
give rise to non-negligible finite-size effects. These ef- 
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FIG. 2. (Color online) Sketch of the spin patterns for several 
of the variational states considered here. The index for the six 
layers of Mn02 in our finite-cluster study is shown at the bot- 
tom. In this figure, only a 2 x 2 unit for each layer is shown. 
The rest of the magnetic state at each layer is obtained by 
repeating the pattern provided. The first four configurations 
(i.e. those with the white background) are uniform and com- 
mon for the n=l and n—2 cases. States of this kind also 
appear in the bulk limit and are not affected by the superlat- 
tice periodicity. The next five configurations shown (with the 
gray background) are only for the n=l case. The last seven 
configurations (with the yellow background) are only for the 
n—2 case. In the non-uniform cases, the cyan bars indicate 
the location of the RO sheets (while the AO sheets are not 
shown explicitly) . The non-uniform phases are denoted using 
a notation mX + m'Y, where X and Y are abbreviations for 
the X-AFM and Y-AFM phases, and m/m' are the fractions 
of the X /Y phases, respectively. 



TABLE I. Parameter profiles used in the present simula- 
tions. Upper rows: the value of Q^iS in the three regions 
of relevance, for the two substrates considered here, namely 
for Mn02 layers in between two RO layers (represented as 
_RMn03), in between a RO layer and a AO layer (represented 
as _R!/2j4i/2Mn03), and finally in between two AO layers 
(represented as AMnOa). The lattice constants of i?Mn03 
and SrTiOs are similar, while those of AMn03 and LaAlOs 
are also similar. Here Q is positive. Lower row: value of Vi 
in the three regions. 

Qs* RMnOs R 1/2 A 1/2 MnQ 3 AMnQ 3 

SrTiOa -Q/2 -Q 

LaAlOa Q Q/2 

Vi V/2 V 



fects are not very prominent for the ferromagnetic (FM), 
G-AFM, or even the CE-type antiferromagnetic cases, 
which were analyzed in previous studies i 25 ' 35 ' 41 However, 
as it will be described below, in the effort described here 
the dominant states involved are the A-AFM and C-AFM 
phases in various combinations. It is necessary to verify 
that the finite-size lattice effects are not too large for 
these phases. Moreover, the presence of a superlattice 
structure may increase the relevance of these effects. 

This present study of size effects and boundary con- 
ditions will be carried out first without the superlattice 
structure, namely mimicking the bulk system. The con- 
clusions of this subsection will then be applied to the 
superlattice study in the rest of the manuscript. In 
Figs. [3] (a-b), the ground state energies (per site) of the 
FM, A-AFM, and C-AFM phases are calculated for the 
case of bulk i? 1/ / 3 A 2 /3Mn03 using the above described 
Hamiltonian (and with Q3i=0 and Vi=0, for simplicity, 
i.e. no superlattice is considered in this test, as already 
explained). The energies of these three states are com- 
pared for both the cases of PBC and TBC. In the PBC 
case, the ground state changes from FM to A-AFM, and 
then to C-AFM with increasing superexchange Jafm (or 
equivalently by decreasing the overall bandwidth). As 
shown in Fig. EJa), with these boundary conditions the 
A-AFM state is stable from Jafm = 0.06 to Jafm = 0.11. 
However, the A-AFM state is actually unstable with the 
use of the TBC, and the ground state changes from FM 
to C-AFM directly, as shown in Fig. \§fb), with the A- 
AFM coexisting with the other two states at just one 
point. Therefore, to reduce this finite-size lattice effect, 
the TBC will be adopted in the rest of the calculations 
described here. However, note that the use of TBC con- 
siderably increases the CPU-time of the study, which ren- 
ders a fully unbiased Monte Carlo simulation too time 
consuming. However, our use of dozens of variational 
states give us confidence that the present results have 
captured the essence of the problem. 

In Figs.FJlJc-d), the double-exchange energies (per site) 
are shown as a function of the value of k in the TBC. 
From the energy values in the vertical axes, it is clear 
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III. RESULTS AND DISCUSSION 
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FIG. 3. (Color online) A comparison of energies using a 4 x 
4x6 lattice (without including the JT distortion Q3; and the 
Coulomb potential Vi, for simplicity) contrasting the cases 
of the PBC and the TBC. The overall e g -electron density is 
chosen as n e =l/3. Panels (a-b) show the total energies (i.e. 
double- exchange and superexchange) per site as a function 
of the superexchange coupling, (a) is for PBC, while (b) is 
for TBC with k=6. The "coincidence" of the three energies 
in (b) at one particular value of Jafm is accidental, (c-e) 
Double-exchange energy for each of the cases indicated, as a 
function of k in the TBC mesh. The accurate energies for the 
bulk limit obtained analytically^ are also shown as dashed 
lines to guide the eye. It is obvious that the use of TBC can 
drastically reduce the finite-size effects. 



that the per-site energy deviation from the bulk value 
using the PBC (corresponding to fc=l) is the most se- 
rious for the case of the C-AFM state (up to — 0.05fo), 
milder for the A-AFM state (up to 0.037irj), and the less 
important for the FM state (difference less than 0.0T9trj). 
And, moreover, fc=6 or even less in the TBC is certainly 
enough for the energy per site to overcome those finite- 
size lattice effects. Thus, this value of k will be adopted 
in the calculations for the superlattices presented below. 



A. (Mn0 3 )„/(AMn0 3 )2„ on SrTiO s 

As explained in the Introduction, the results of our 
investigations will be provided for two particular sub- 
strates, starting with SrTiOs. For this case, note that 
the average lattice constant for a pseudo-cubic i?MnC>3 
(i?=La, Pr, Nd) is close to the lattice constant of SrTiOs 
(3.905 A). More specifically, this (average) lattice con- 
stant is about 3.94 A for i?=La, 3.916 A for i?=Pr, and 
3.907 A for i?=Nd. 43 Then, in the numerical simulations 
to describe this particular case, the ilMnOa layers grown 
on SrTiOa are assumed to be nearly cubic, namely the 
value Qsi = is used for the active JT mode in our study. 
In contrast, the cubic AMnO^ has a much smaller lattice 
constant, in particular 3.805 A for A=Sr and 3.727 A 
for A=C&r^ which gives rise to an important in-plane 
expansion and out-of-plane reduction of the lattice con- 
stants when growing on the SrTiOs substrate. The ratio 
c/a of AM11O3 on SrTiOs will be smaller than 1, implying 
a negative Q^i- Therefore, in the following simulations 
with SrTiOa as the substrate, the values of the Q^ vari- 
ables are considered to be in the i?Mn03 region (more 
specifically, when a Mn02 layer is in between two RO 
layers), but — Q in the AVhiOz region (i.e. when the 
Mn02 layer is in between two AO layers), where Q is a 
positive parameter to be tuned in our study. For the in- 
terfacial i?0-Mn02-AO layer, an average value —Q/2 is 
used as a first order approximation, as shown in Table HI 

The DFT results^SS on (LaMn0 3 )2 n /(SrMn0 3 )„ and 
(LaMn03)„/(SrMn03)2n superlattices (n=l,2) have in- 
dicated that the Coulomb energy difference between 
LaMn03 and SrMnOs is about 0.7-1 eV, which cor- 
responds to 1.4-2^0 in our simulation if to is 0.5 eV 
(a reasonable approximate value for wide bandwidth 
manganites). 25 However, besides the Coulomb potential 
from the A-site cations, the V, used here also contains 
other contributions, such as those from the breathing 
phononic Q\ mode. In fact, previous model studies found 
that T^=0.6-0.9to is a proper range for this potential in 
the double-exchange model. 25 However, for completeness 
here the full range of V up to 2irj will be investigated in 
the phase diagrams. 

The ground state phase diagrams corresponding to 
(i?Mn0 3 )i/(AMn0 3 ) 2 (n=l) using SrTi0 3 as the sub- 
strate are shown in Fig. [4] at various values of XQ. With- 
out the JT distortion (i.e. for XQ — 0, see Fig. @] (a)), 
the A-AFM phase is not stable in this superlattice despite 
varying the parameters Jafm and V. Instead, FM ten- 
dencies dominate with at most just 1/3 of the superlat- 
tice being G-AFM even at the largest Jafm studied here. 
The presence of one G-AFM layer for every two FM lay- 
ers is in agreement with the naive intuition described in 
the Introduction, namely that the G-AFM state should 
be favored inside AMn03. However, these results change 
qualitatively when the strain effects are incorporated. In 
fact, with increasing XQ, the A-AFM phase now emerges 
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FIG. 4. (Color online) Phase diagram for the case of 
(_RMn0 3 )i/(AMn0 3 )2 with SrTi0 3 as substrate, at several 
values of XQs. 



FIG. 5. (Color online) Phase diagram for the case of 
(i?Mn0 3 ) 2 /CAMn0 3 ) 4 with SrTi0 3 as substrate, at several 
values of XQ. 



around Jafm=0.05-0.06£o and takes over an increasing 
fraction of parameter space with increasing XQ. This ef- 
fect is the most prominent in the absence of the electro- 
static potential. In principle, a strong Coulomb potential 
V, which reduces the e g electronic density in the AM11O3 
region, can suppress this uniform A- AFM order. In fact, 
in Fig. |4] (b) the A- AFM phase is suppressed rapidly with 
increasing V. However, in spite of this negative influ- 
ence of V, when the JT distortion is strong enough, e.g. 
XQ ^ OMo, the A-AFM phase is stable and quite robust 
in the broad region of parameter space investigated here, 
in agreement with the experiments described in the Intro- 
duction. Note also the appearance of additional robust 
states with further increasing Jafm, particularly the Cx 
state. As discussed later in the manuscript, these states 
could appear in manganites with a bandwidth smaller 
than for the case A=Sr. 

For the case of the (i?Mn03)2/(^4Mn03)4 superlattice 
with SrTiC>3 as substrate, the strong JT distortion is 
even more crucial to obtain the A-AFM phase, as shown 
in Fig. [5] The A-AFM phase is stable only in the interval 
Jafm=0. 05-0. 06^o for the cases of a strong JT distortion 
and weak to intermediate Coulomb potentials. Compar- 
ing with the previously described case of n = 1, the A- 
AFM phase appears to be more fragile in the n = 2 case, 
with stability regions in parameter space that are smaller 
in size for n=2 than for n=l. This tendency is in agree- 
ment with the experimental results reporting that the 
A-AFM Neel temperature of (LaMn0 3 )i/(SrMn0 3 ) 2 is 
higher than (LaMn0 3 )2/(SrMn0 3 )4. 21 The physical rea- 
son is that the charge density deep into the AMn03 re- 



gion is very low due to the Coulomb potential V . Thus, 
the combination of 3 FM layers (i?Mn03 and interfacial 
layers) plus 3 G-AFM layers {AM11O3 layers) eventually 
wins in the energy competition when V is large enough, 
in agreement with the naive expectations presented in the 
Introduction. However, at intermediate values of Jafm 
and for large enough strain, represented by XQ, the ex- 
perimentally observed A-AFM state is stabilized. 

Increasing further Jafm, the ground state becomes a 
Cx-type antiferromagnetic (Cx-AFM) phase. Here, the 
Cx-AFM is similar to the C-AFM phase, but the chains 
with the spins aligned lay in the a — b plane while in the 
standard C-AFM state this spin alignment is along the 
c axis (out-of-plane). This phase is certainly interesting 
but our description of results here and in the rest of the 
publication will focus mainly on the phases that appear 
at intermediate values of Jafm- The reason is that for 
the case of bulk manganites, phases not considered in 
our effort such as the CE-AFM state with Q2 distortions 
also appear in the range of Jafm studied here at least at 
particular electronic densities such as half-doping. More- 
over, other states stabilized at larger values of Jafm arid 
the same doping, such as the G-AFM state, have not 
been found experimentally for bulk compounds. Thus, it 
is natural that the focus of our results for superlattices 
should only be on the intermediate Jafm range. Never- 
theless, a brief discussion for narrow bandwidth mangan- 
ites is presented below in this manuscript as well. 
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FIG. 6. (Color online) Phase diagram for the case of 
(7?Mn03)i/(ylMn03)2 with LaA103 as substrate, at several 
values of XQ. 



FIG. 7. (Color online) Phase diagram for the case of 
(_RMn03)2/(^4Mn03)4 with LaA103 as substrate, at several 
values of XQ. 



B. (EMn0 3 )„/( J 4Mn0 3 )2„ on LaA10 3 

Consider now the case of LaAlC>3 as the sub- 
strate. In a recent first-principles calculation, the 
(LaMn03) n / (SrMn03)2n superlattices grown on LaAlC>3 
as the substrate were also studied, and results quite dif- 
ferent from those of the SrTiC>3 substrate were obtained, 
as already mentioned in the Introduction^ For this rea- 
son, here the phase diagrams for superlattices grown on 
LaA103 are also calculated to contrast model Hamilto- 
nian results against first-principles approaches. 

In contrast to the case of SrTiC>3 as substrate, the 
lattice constant of LaAlC>3 (which is approximately 3.791 
A) is much smaller than the lattice constant of -RM11O3, 
but it is close to SrMn0 3 (3.805 A) and to CaMn0 3 
(3.727 A). Therefore, the i?MnC>3 layers are considerably 
compressed in-plane while the AM11O3 layers are almost 
free from strain when LaA103 is the substrate. Then, 
to simulate this strain effect a positive value of Q^i = 
Q is adopted for the i?MnC>3 layers while Qzi = is 
used in the ^4Mn03 region. Similarly as in the previous 
subsection, for the interfacial layers Qzi is taken as Q/2 
as a first order approximation, as shown in Table HI This 
simple and smooth JT profile is the opposite to that used 
before for the SrTiC-3 case. 

As shown in Fig. [51 the A- AFM phase, which is ro- 
bust in the phase diagrams for the case of SrTiOs as 
substrate, no longer survives for the case of LaAlC>3. In- 
stead, now the C-AFM state becomes the most likely 
phase appearing in reasonable regions of parameter space 
for the large bandwidth manganites. This result agrees 
with the first-principles prediction that the C-AFM phase 



is quite robust in the (LaMn03)/(SrMn03)2 superlat- 
tice on LaA103^ Note that in Fig. [5] strain is crucial to 
stabilize the C-AFM phase. Without strain, the naive 
expectation described in the Introduction of finding G- 
AFM layers in the ^4Mn03 region is indeed fulfill since 
the 2FM+1G and 2FM+1G:B states dominate. Once 
again, the crucial role of strain in these investigations 
becomes clear. 

Note, however, that for (i?Mn0 3 ) 2 /(^Mn03) 4 the uni- 
form C-AFM phase appears to be fragile, and becomes 
stable only in a small region varying V (Fig. [7|). In- 
stead, the hybrid configuration consisting of 3 FM lay- 
ers (i?MnC>3 and interfacial layers) plus 3 C-AFM layers 
(AMn03 layers) is mainly stabilized when V increases. 
Thus, the experimental realization of the C-AFM state 
for n=2 using LAO as substrate appears more problem- 
atic than for n=l. Further increasing n, eventually the 
naive expectation of a G-AFM dominance in the AMnOs 
region must be satisfied, but it is difficult to predict the 
precise value of n for the transition from a C-AFM dom- 
inated to a G-AFM dominated behavior in the AM11O3 
layers. 



C. Electronic density and orbital occupation 

The e g electronic density and the orbital occupan- 
cies are important factors to understand intuitively the 
origin of the previously described ground states of the 
(i?Mn03)„/(^4Mn03)2n superlattices on different sub- 
strates. These quantities are closely related to both the 
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FIG. 8. (Color online) The e g density profiles of the orbital 
levels in (_RMn03)n/(j4Mn03)2n for the special example of 
l / =0.6fo and AQ=0.9to. In all panels, the values of n and 
the substrate used are indicated. The total e g density (black) 
and the orbital occupations (blue: d x ?_ y i\ red: d iz i_ r i) are 
shown. In (c) and (d), the d- iz i_ r i orbital is virtually the 
only occupied orbital, and for this reason the total density is 
not shown since it is virtually identical to the d 3z 2_ r 2 density. 
The cyan bars denote the RO sheets in the superlattices. 
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FIG. 9. (Color online) Sketch of the dominant orbital occu- 
pation for the two substrates and the two values of n studied 
in this manuscript. The six drawings for each case correspond 
to the six layers studied here, and their size is proportional 
to the actual occupation of the respective orbitals. The lines 
in cyan color indicate where the RO layers are. Since the 
associated A-AFM and C-AFM magnetic orders are uniform 
in the layers that form the superlattice, then the results in 
this figure are the same for all the other chains running in the 
direction perpendicular to the interfaces. 



uniform Q3 JT distortions (determined by the substrates) 
and the spin arrangements. 

As a typical example, the case of V~0.6to and 
AQ=0.9£o will be discussed in detail, and the e g electronic 
densities and their orbital components will be calculated 
explicitly. As shown in Figs. [5] (a) and (b), for the case of 
the A-AFM phase on SrTiC>3 the e g densities only weakly 
fluctuate around the average value n e =l/3 all through 
the six layers of the superlattice studied here, despite 
the presence of the electrostatic potential modulation. 
In fact, even the lowest e g density within the AM11O3 
region is still quite robust and only slightly smaller than 
0.3, i.e. it is appreciably higher than the e g electronic 
density of AMnC>3 in the bulk (i.e. ~ 0). In other words, 
the transfer of charge from i?Mn03 to AMn03 in these 
superlattices with thin components is quite prominent 
and important to stabilize the particular spin orders that 
have been discussed in this publication. Therefore, the 
magnetic orders found in the present manganite super- 
lattices cannot be simply extended to pure i?Mn03 and 
AMnC>3 films even under the same strain conditions. 

This is a crucial observation to understand the sta- 
bilization of the A-AFM order in the superlattice 
as compared to the G-AFM order in the bulk of 
AMn03. Note that in previous theoretical studies of 



(LaMn0 3 )2n/(SrMn0 3 ) n superlattices^^ the G-AFM 
order was found to be stable in the SrMn03 region be- 
cause its local electron density was much lower than 
found in the example studied here. The reason for this 
difference is the inclusion of strain effects (i.e. a nonzero 
XQ) in the present analysis. In fact, the effort described 
here suggests that it may be necessary to revisit the 
previously studied (LaMn03)2n/(SrMn03)„ system, be- 
cause the effects of strain were not taken into account in 
those previous theoretical investigations. 

For the case of the C-AFM phases on LaA103 the sit- 
uation is different. In this case, the orbital occupation 
is purely d 3z 2_ r 2 1 without an appreciable d x i_ y i compo- 
nent. The e g electronic densities are also substantially 
modulated, following much closer the cations than for 
the case of SrTi03. This charge modulation is prominent 
for the C-AFM phases on LaA10 3 , while for the A-AFM 
phases on SrTiOs the electronic density is almost uni- 
form for n—1 and only weakly modulated for n—2, as it 
is clear by comparing Figs. [8] (a)-(b) with (c)-(d). 

Addressing now the orbital occupation, for the A-AFM 
phase on SrTi0 3 the d x i_ y i is the orbital the most popu- 
lated, with a mild contribution of the d^ z 2_ r i component 
in the i?Mn03 portion of the superlattice for n—2, as 
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FIG. 10. (Color online) The e g -orbitals energy splitting and 
the electrostatic potentials in (i?Mn03)„/(ylMn03)2n for the 
special example of V—0.6to and AQ=0.9to- In all panels, the 
values of n and the substrate used are indicated. The final 
energy levels of the two orbitals (after the JT splitting plus 
the shift from Vi) are shown as the blue (d x 2_ y i) and red 
(d 3z 2_ r 2 ) solid curves. The cyan bars denote the RO sheets in 
the superlattices. The electrostatic potentials Vi are shown in 
dashed black curves. For completeness, the energy split of the 
two orbitals only due to the JT Q3 mode is also shown, with 
the blue (red) dashed curve denoting the d x 2_ y 2 (d 3z 2_ r 2) 
orbital. 

shown in Figs. [HI (a, b). On the other hand, for the C-AFM 
phase on LaA103, it is the d^ z 2_ r 2 orbital that dominates, 
as intuitively expected (see Figs. [9] (c,d)). Figs. [9] (a-d) 
show only one line of orbitals running perpendicularly to 
the interfaces because the orbital occupation is the same 
in all the sites of the x-y layers since the A-AFM and 
C-AFM are uniform in those layers. 

To understand better the density and orbital profiles 
found in our variational study, the splitting of the energy 
between the two e g -orbitals induced by the Q3 mode is 
shown in Figs. [TU] (a-d). On SrTiOa, the d x 2_ y 2 and 
d 3z 2_ r 2 orbital levels in the AMn03 region of the su- 
perlattice are shifted down and up, respectively, due to 
the influence of the XQ^i coupling in that region, while 
they are degenerate for n=2 in the i?Mn03 portion of 
the superlattice due to its nearly cubic lattice symme- 
try (\Q 3i =0 there). Note that for n=l, XQ 3i is always 
nonzero, thus there is always a splitting between the two 
orbitals. Therefore, the d x 2_ y 2 level in the AMn03 re- 
gion is low enough to "accumulate" e g electrons even 
when the electrostatic potential Vi is imposed, as shown 



in Figs. [10] (a,b). Since the d x 2_ y 2 orbital, with its in- 
plane lobes, prefers the in-plane double-exchange hop- 
ping between NN sites, then the A-AFM state is sta- 
bilized since in this case the double-exchange process is 
precisely restricted to occur within the x-y plane. 

Contrary to the case of SrTi03, when LaA10 3 is the 
substrate the d x 2_ y 2 and d 3z 2_ r 2 orbital levels in the 
i?Mn03 portion of the superlattice are shifted up and 
down, respectively, due to the influence of XQsi, while 
they are degenerate in the AMn03 component of the 
superlattice where XQ3i—0. Since the double-exchange 
process in the C-AFM phase is only active along the out- 
of-plane chains with the same orientation of the spin and 
since in this direction the hoppings associated with the 
x 2 — y 2 vanish (t z aa =tl a =t z ab =Q) , then the d 3z 2_ r 2 orbital 
is occupied while the d x 2_ y 2 orbital is virtually empty. 

In summary, the two e g orbitals are both active in un- 
strained layers but typically only one dominates in the 
uniformly strained layers. This simple observation is suf- 
ficient to understand several features of the phase dia- 
grams that were obtained numerically. 



D. Extensions to narrow-bandwidth manganite 
superlattices 

Most experiments on Mn-oxide superlattices have fo- 
cused on the wide bandwidth manganites involving 
LaMn0 3 and SrMn0 3 as components Mr-^^- For this 
reason considerable theoretical investigations have been 
carried out with this focus as wel h 25 i 26 ' 28 However, 
some recent theoretical efforts pointed out that exotic 
phenomena might exist at interfaces involving narrow- 
bandwidth manganites, even including the presence of 
novel magnetic states that do not appear in the bulk 
phase diagrams of the individual materials that form the 
superlattice. 35 Since reducing the bandwidth is equiva- 
lent to increasing the values of Jafm and V, because 
they are in units of the hopping to, then the previously 
discussed phase diagrams contained in Figs. 0][7] indicate 
that other magnetic orders beyond the uniform A-AFM 
and predicted C-AFM states could be stabilized in these 
superlattices if the parameters Jafm or V are increased 
further by using manganites different from LaMnOs or 
SrMn0 3 . 

For example, considering the case of the SrTiOs sub- 
strate, the uniform Cx-AFM order is neighboring the A- 
AFM order in the phase diagrams Figs. 0] and [SJ and it 
appears to be quite robust in both the n = 1 and n = 2 
superlattices for a sufficiently large Jafm (correspond- 
ing to narrow bandwidth cases). It should be noted that 
the Cx-AFM order does not exist in bulk form, thus its 
stabilization in superlattices would be quite interesting. 
The Cx-AFM phase, which breaks the in-plane symme- 
try between the two directions, will result in macroscopic 
physical anisotropies in, e.g., transport properties even 
when the underlying substrate is cubic. 

The phase diagrams for the case of SrTiOs as sub- 
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strate include several other phases besides the A-AFM 
and Cx-AFM states. These other phases correspond to 
combinations of FM, A-AFM, and G-AFM states, with 
different patterns, as shown in Figs. 2] and [5] For config- 
urations involving a FM component, note that it is quite 
difficult to obtain a macroscopic magnetization in these 
superlattices, since there is no place for FM order in the 
phase diagrams for manganites when in bulk form and 
for hole doping x > O.5. 27 ' 37 Thus, if a net magnetic mo- 
ment could be stabilized via superlattices for a nominal 
composition x larger than 0.5, this would also represent 
a new state that is not present in bulk form. Moreover, 
combinations of FM and G-AFM layers could potentially 
be used as spin valves by controlling the local magnetic 
orders with magnetic fields. 

With regards to the LaA10 3 substrate, the C-AFM 
state is quite robust in the n=l case (Fig. [6]), which may 
imply an enhanced Neel temperature, analogous to the 
enhanced A-AFM Neel temperature observed experimen- 
tally for SrTi0 3 as substrate. However, in the n=2 case 
the C-AFM state occupies a smaller portion of the phase 
diagram and states such as those involving 3 FM layers 
plus 3 G-AFM layers may become stable if V is increased 
further, similarly as in the case using SrTi03. Other 
combinations such as the 5 C-AFM layers plus 1 G-AFM 
layer are also possible if the bandwidth is reduced. In 
this case, and also for SrTi03, the phase diagrams are 
rich, implying states that are very close in energy, and it 
is possible that external magnetic fields could be used to 
tune the ground state properties of these superlattices. 

IV. SUMMARY 

Summarizing, here the two-orbital double-exchange 
model was employed to investigate the phase diagrams 



of the (i?Mn0 3 ) n /(^4Mn0 3 ) 2 „ superlattices. Two differ- 
ent substrates have been tested to illustrate two different 
types of strain effects, giving rise to distinct ground states 
and phase diagrams. The underlying dominant physical 
mechanism is the uniform splitting of the e g orbital lev- 
els by the Q 3 mode of Jahn- Teller distortion imposed by 
the substrate strain. Besides the uniform A-type and C- 
type antiferromagnetic phases that are stabilized in the 
intermediate range of couplings for the case of the wide- 
bandwidth manganite superlattices, our study also re- 
veals several other possible phases, including many with 
non-uniform modulations of the spin order. These ex- 
otic patterns could be found in superlattices built using 
narrow-bandwidth manganites, an area of research that 
should be more actively pursued experimentally due to 
its potential to unveil new states of manganites. 
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